Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file vl.4) 



The Lx — T Relation and Intracluster Gas Fractions of 
X— ray Clusters 



ON 

c 

on 



Monique Arnaud 1 and August E. Evrard 2 ' 3 

1 C.E.A., DSM, DAPNIA, Service d'Astrophysique, C.E. Saclay, F-91191, Gif-Sur-Yvette Cedex, France 

2 Institut d'Astrophysique, 98bis Blvd Arago, 75014 Paris, France 

3 Physics Department, University of Michigan, Ann Arbor, MI 48109-1120 USA 
arnaud@sapvxg.saclay.cea.fr 

evrard@umich. edu 



Accepted 1999 January 12 



(N 
> 
CO 

in 

CO 
NO 
O 
OO 
On 

^ 

Of 

6 



ABSTRACT 

We re-examine the X-ray luminosity— temperature relation using a nearly homo- 
geneous data set of 24 clusters selected for statistically accurate temperature measure- 
ments and absence of strong cooling flows. The data exhibit a remarkably tight power- 
law relation between bolometric luminosity and temperature with a slope 2.88 ± 0.15. 
With reasonable assumptions regarding cluster structure, we infer an upper limit 
on fractional variations in the intracluster gas fraction ((Sf gas / fgas) 2 ) 1 ^ 2 < 15%. A 
strictly homogeneous GINGA subset of 18 clusters places a more stringent limit of 
9%. 

Imaging data from the literature are employed to determine absolute values of f gas 
within spheres encompassing density contrasts S c = 500 and 200 with respect to the 
critical density. Comparing binding mass estimates based on the virial theorem (VT) 
and the hydrostatic, /3-model (BM), we find a temperature-dependent discrepancy 
in fgas between the two methods caused by systematic variation of the outer slope 
parameter (3 with temperature. Mean values (for H o = 50 km s" 1 Mpc" 1 ) range from 
f gas = 0.10 for cool (T < 4 keV) clusters using the VT at 5 C = 500 to 0.22 for hot 
(T > 4 keV) clusters using the BM at 5 C = 200. There is evidence that cool clusters 
have a lower mean gas fraction than hot clusters, but it is not possible to assess the 
statistical significance of this effect in the present dataset. The T dependence of the 
ICM density structure, coupled with the increase of the gas fraction with T in the VT 
approach, explains the steepening of the Lx~T relation. 

The small variation about the mean gas fraction within this majority sub- 
population of clusters presents an important constraint for theories of galaxy formation 
and supports arguments against an Einstein-deSitter universe based on the population 
mean gas fraction and conventional, primordial nucleosynthesis. The apparent trend 
of lower gas fractions and more extended atmospheres in low temperature systems are 
consistent with expectations of models incorporating the effects of galactic winds on 
the intracluster medium. 
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1 INTRODUCTION 

The relation between X-ray temperature T and luminos- 
ity Lx is a sensitive diagnostic of structural regularity in 
clusters of galaxies. Because the luminosity is governed by 
the mass of gas in the intracluster medium (ICM) while the 



temperature is determined by the total, gravitating cluster 
mass, the Lx~T relation can be used as a tool to probe vari- 
ations in gas fraction — the ratio of ICM gas to total mass 
— in these systems. 

Evidence for a correlation between these basic observa- 
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tional quantities has existed since the early days of X-ray 
astronomy (Mitchell et al. 1977, 1979; Mushotzky et al. 1978; 
Edge & Stewart 1991, David et al. 1993). Despite its rela- 
tive maturity, a key characteristic of this relation was only 
recently established. Fabian et al. (1994) demonstrate that 
the departure of a cluster from the mean Lx — T relation 
is correlated with the strength of the emission associated 
with its cooling flow core. The picture resulting from this 
work is of a mixed population which is likely to be time 
variable. Fabian et al. speculate that cooling flows are a re- 
current phenomenon associated with a secular instability in 
the ICM plasma which periodically is interrupted and re- 
set by strong merging encounters. Periodic mergers of vary- 
ing strength are a natural feature of hierarchical clustering 
models of structure formation, and numerical simulations 
suggest that cooling flow features can be erased by merg- 
ers (Garasi, Burns & Loken 1997). However, the theoretical 
picture remains incomplete because of the complexities of 
the physics operating in the cluster core. Regardless of the 
exact physical mechanisms responsible, the empirical fact is 
that the cluster population can be broadly classed into two 
categories by making a cut in cooling flow strength M. 

In this paper, we present and analyze the luminosity- 
temperature relation from a nearly homogeneous data set 
of weak cooling flow clusters, where M < 100 Mq yr _1 de- 
fines weak. The data set is described in §2, with particular 
attention paid to systematic errors in determination of both 
Lx and T. In §3, we present the Lx~T relation and show 
that these data place stringent limits on gas fraction varia- 
tions within this sub-population. Imaging data taken from 
the literature allow us to estimate absolute values of the gas 
fraction at fixed density contrast in §4. We compare results 
for two different binding mass estimation methods, and em- 
phasize the systematic uncertainty introduced by this choice. 

Hubble constant dependencies are displayed in the pa- 
per via /150 = H o /50 km s _1 Mpc - and we assume h^ — l 
when quoting numerical values throughout the paper. 



2 THE SAMPLE 

We analyse an archival data set, part of which was assembled 
and discussed by Arnaud (1994). A prime objective is to 
limit both statistical and systematic errors in L x and T, 
in order to obtain an accurate assessment of the intrinsic 
dispersion about the mean L x — T relation. 



2.1 Selection Criteria 

Inclusion of a cluster in the sample is based on three cri- 
teria : (i) small statistical errors in measured temperature 
AT/T<10% (at 90% confidence), (ii) weak or absent cool- 
ing flow M< 100 M yr -1 , and (iii) kT > 2 keV. The value 
of M is set by the relative importance of core luminosity to 
the total, and is typically evaluated using models of multi- 
phase, steady-state accretion (Thomas, Fabian & Nulsen 
1987). Since nearly two-thirds (33/51) of the clusters in the 
X-ray flux limited sample of Edge et al. (1990) have flows 



below this limit (Fabian et al. 1994), our analysis is of a 
majority population. 

Table 1 presents a listing of the 24 clusters satisfying the 
above criteria. Eighteen have both T and Lx determined by 
the GINGA satellite (references are given in the table). Lu- 
minosities for three nearby clusters (Virgo, in particular) 
are calculated from scanning data, which accounts directly 
for the extended emission outside the 1.1 x 2 sq deg FOV 
(FWHM) of the satellite instrument. Temperature measure- 
ments for the other six clusters come from ASCA, with lu- 
minosities from Einstein MPC (3), EXOSAT (2), and ASCA 
(!)• 

We believe this list to be complete with respect to publi- 
cations available end of 1996. Sixteen clusters in our sample 
are included in the X-ray flux limited sample of Edge et al. 
(1990). Two more clusters (A2163 and A1413) have fluxes 
above the Edge et al flux limit (1.7 x 10~ n erg s _1 cm -2 ). 
Our sample is ~ 70% complete as compared to the Edge 
et al. sample down to a flux of 3.4 x 10~ n erg s _1 cm~ 2 
(considering from now on only weak cooling flow clusters 
in the Edge sample for consistency) and highly incomplete 
below this flux (3 clusters among 20). The other clusters of 
our sample below the Edge et al. flux limit are either distant 
luminous clusters (A3921, A665, A2218, A370) or poor clus- 
ters (A2634 and AWM4). The temperature distribution of 
our sample shows a clear deficit of clusters at intermediate 
temperature (4 — 6 keV) as compared to the temperature 
distribution of the Edge et al. sample. 

Included in Table 1 are the measured temperature and 
luminosity in the 2-10 keV energy band (cluster rest frame 
and corrected for absorption). We convert the 2-10 keV 
luminosity to bolometric luminosities using an isothermal 
plasma emission model (Mewe, Gronenshild & van den Oord 
1985; Mewe, Lemen & van den Oord 1986) and the measured 
temperatures and abundances. The bolometric correction 
factor depends mostly on the temperature; it is higher at 
low temperature (~ 3.4 at 2.2 keV), decreases with increas- 
ing temperature up to about 8 keV (~ 2.1) and increases 
again above (~ 2.4 at 15 keV). The statistical errors on the 
luminosities are always much smaller than the correspond- 
ing errors on the temperature and are thus not taken into 
account here. 

2.2 Systematic Errors 

To estimate accurately the intrinsic dispersion of the Lx — T 
relation requires an attempt to understand and quantify 
sources of systematic errors. Such errors can arise from 
calibration uncertainties, systematic errors in background 
subtraction and, for luminosities deduced from collimated 
instruments (EXOSAT, GINGA), improper correction for 
cluster extent (loss of efficiency at large radii, emission out- 
side the FOV) . Further uncertainties due to the uncertainties 
on the plasma emission model are negligible for clusters in 
the considered temperature range (see Arnaud et al. 1991; 
Arnaud et al. 1992). 

We attempt to guage the magnitude of systematic un- 
certainty in the bolometric luminosity by comparing to val- 
ues derived from the ROSAT All-Sky Survey (RASS) by 
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Table 1. Basic data for the 24 clusters in the sample 



Name Rcdshift Lx [2-10 kov] Temperature Abund. Lbolom. P Core radius Ref. 

(ergs/s) (kcV) (rel. solar) (crgs/s) (kpc) 



Virgo 


0.0038 


2.00 


10 43 




2.20 ±0.02 


0.44 


6.81 


10 43 


0.46 


12.6 


29,29,7,7 


A262 


0.0164 


3.70 


10 43 




2.41 ±0.05 


0.55 


1.15 


10 44 


0.53 ± 0.03 


95 ± 25 


34,34,12,19 


A644 


0.0704 


1.06 


10 45 




6.59 ±0.17 


0.33 


2.24 


10 45 


0.70 


204 


16,16,20,20 


A754 


0.0534 


1.16 


10 45 




7.57 ±0.31 


0.23 


2.46 


10 45 






1,1 


A1367 


0.0215 


9.09 


10 43 




3.50 ± 0.18 


0.31 


2.26 


10 44 


0.50 ± 0.10 


390+? 1 ? 


16,16,29,20 


A1656 


0.0232 


1.13 


1Q 45 




8.21 ±0.09 


0.22 


2.41 


10 45 


0.75 ± 0.03 


420 ± 24 


18,18,8,8 


A2256 


0.0601 


1.11 


10 45 




7.51 ±0.19 


0.28 


2.35 


10 45 


0.81 ± 0.01 


542 ± 12 


16,16,9,9 


A2319 


0.0564 


2.24 


10 45 




9.12 ±0.15 


0.24 


4.83 


10 45 


0.68 ± 0.05 


410 ± 50 


34,34,19,19 


A2634 


0.0321 


7.09 


10 43 




3.07 ±0.11 


0.38 


1.89 


10 44 


0.58 


320 


1,1,11,11 


A3526 


0.0109 


7.00 


10 43 




3.54 ±0.13 


0.54 


1.70 


10 44 


0.45 ± 0.03 


93 ± 19 


32,33,23,23 


A3921 


0.0940 


5.20 


10 44 




4.90 ±0.55 


0.26 


1.15 


10 45 


0.75 ± 0.04 


385 ± 42 


2,2,2,2 


AWM4 


0.0315 


2.53 


10 43 




9 4^ -1- 98 


0.93 


7.52 


10 43 


0.43 


190 


31 31 31 31 


AWM7 


0.0179 


1.70 


10 44 




3.82 ±0.15 


0.49 


3.99 


10 44 


0.53 ±0.01 


102 ±5 


31,31,28,28 


Ophui. 


0.0280 


1.90 


10 45 




9.10 ±0.30 


0.24 


4.10 


10 45 


62+ ' 01 


177t« 


21,21,24,24 


A665 


0.1820 


1.50 


10 45 




8.26 ±0.90 


0.49 


3.16 


10 45 


0.66 ± 0.05 


379 ± 70 


17,17,5,5 


A1413 


0.1427 


1.82 


10 45 




8.85 ±0.50 


0.19 


3.92 


10 45 


0.62 


156 


16,16,11,11 


A2163 


0.2010 


6.00 


10 45 




14.6 ± 0.85 


0.40 


1.45 


1Q 46 


0.62 ± 0.02 


305 ± 19 


14,14,14,14 


A2218 


0.1710 


9.40 


10 44 




6.70 ±0.50 


0.20 


2.00 


10 45 


fi5+ - 08 
u - oo -0.05 


2261*° 


25,25,6,6 


A370 


0.3700 


1.30 


10 45 


[1] 


8.80 ±0.80 


0.50 


2.76 


10 45 






4,4 


A399 


0.0715 


7.76 


10 44 


[3] 


7.40 ± 0.50 


0.25 


1.64 


10 45 


0.52 ± 0.05 


215 ±35 


13,15,19,19 


A401 


0.0748 


1.47 


10 45 


[3] 


8.40 ±0.50 


0.21 


3.15 


10 45 


0.61 ± 0.01 


285 ± 10 


13,15,19,19 


A1060 


0.0114 


2.50 


10 43 


[2] 


3.10 ±0.20 


0.40 


6.60 


10 43 


0.61 


94. 


13,30,22,22 


A3558 


0.0478 


4.23 


10 44 


[2] 


5.50 ±0.25 




9.06 


10 44 


0.61 


324 


13,27,3,3 


Trian. 


0.0510 


1.26 


10 45 


[2] 


10.3 ±0.80 


0.26 


2.77 


10 45 


0.63 ± 0.02 


286 ± 16 


13,26,26,26 



Notes: column (9) references for the data (respectively X-ray luminosity, temperature , j3 and core radius): 1. Arnaud 
(1994); 2. Arnaud ct al. (1997); 3. Bardclli et al. (1996); 4. Bautz ct al. (1994); 5. Birkinshaw, Hughes & Arnaud (1991); 6. 
Birkinshaw & Hughes (1994); 7. Bohringer et al. (1994); 8. Briel, Henry & Bohringer (1992); 9. Bricl & Henry (1994); 10. 
Buotc & Canizarcs (1996); 11. Cirimele, Nesci & Trevese (1997); 12. David, Jones & Forman (1995); 13. Edge et al. (1990); 
14. Elbaz, Arnaud & Bohringer (1995); 15. Fujita et al. (1996); 16. Hatsukade (1989); 17. Hughes & Tanaka (1992); 18. 
Hughes et al. (1993); 19. Jones & Forman (1984); 20. Jones & Forman (1990, priv. comm.); 21. Kafuku ct al. (1992); 22. 
Loewenstein & Mushotzky (1996); 23. Matilsky, Jones & Forman (1985); 24. Matsuzawa et al. (1996); 25. McHardy et al. 
(1990); 26. Markevitch, Sarazin & Irwin (1996); 27. Markcvitch & Vikhlinin (1997); 28. Neumann & Bohringer (1995); 29. 
Takano (1990) and Koyama (1997, priv. comm.); 30. Tamara ct al. (1996); 31. Tsuru (1993); 32. Yamanaka & Fukazawa 
(priv. comm. in Fukazawa et al. (1994)); 33. Yamashita 1992; 34. Yamashita (1997, priv. comm.) 



Ebeling et al. (1996) for Abell Clusters. Eigthteen of the 
clusters in our sample are in this dataset. We convert the 
RASS band-limited, total unabsorbed flux, derived from a 
procedure described in the above reference, to bolometric 
Lx using the same plasma emission model as that used on 
our own data. In addition, we include RASS data for Virgo 
from Bohringer et al. (1994) and ROSAT imaging data for 
Triangulum (Markevitch, Sarazin & Irwin, 1996) and Ophi- 
uchus (Buote & Tsai 1996), from which we calculate total 
luminosities by extrapolating the /3-model fits given in the 
references. Conversion to a bolometric measure for imaging 
data is the same as for the RASS data. 

We find good agreement between the GINGA and 
ROSAT bolometric luminosities; 16 objects yield a mean ra- 
tio of 1.00 with standard deviation of 0.19. The agreement is 
not as good for the five non-GINGA members of our sample, 
but the difference is driven entirely by two clusters — Tri- 
angulum and A3558 — for which the ROSAT luminosities 



are larger than the EXOSAT values we use by factors of 1.72 
and 1.56, respectively. Both clusters are nearby and have ex- 
tent comparable or larger than the EXOSAT FOV (45' x 45' 
FWHM), the virial radius (defined below) of these 2 clusters 
are respectively at 49 and 38 arcmin from the center. Their 
flux can be underestimated due to the neglect of the loss 
of efficiency at large radii. Furthermore, both observations 
were significantly offset, the pointing offset being specially 
important for A3558 (38 arcmin, Edge 1989). Although the 
fluxes published by Edge et al. (1990) include correction for 
this effect (correction as for a point source), this is likely 
to induce further uncertainties in the flux estimate. For the 
sake of consistency with the remaining dataset, we applied 
no correction to the published values. 

The comparison of the GINGA and ROSAT luminosi- 
ties, which constitutes the majority of our sample, indicates 
that the scatter expected from systematic uncertainty in Lx 
is not larger than 20% (rms). Clusters with large angular 
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extent — in particular, larger than GINGA's 1.1 deg x2deg 
field of view — are likely to incur larger error, but the num- 
ber of such objects in our sample is small. Furthermore, we 
looked for systematic variation with redshift of the ratio be- 
tween the GINGA and ROSAT luminosities. Such a depen- 
dence with redshift was noticed by Ebeling (1993) for the 
ratio between the luminosities (HE AO, MPC and EXOSAT 
data) given by Henry & Arnaud (1991) and the RASS lu- 
minosities. The observed increase of this ratio with redshift 
(up to z ~ 0.07) was interpreted as being due to an un- 
derestimate of flux for nearby, extended clusters by earlier 
experiments. This effect is not apparent in the GINGA data, 
due to the larger instrument FOV and the use of scanning 
mode for the most extended clusters. 

We now consider systematic errors on temperature. A 
detailed estimate of this error, due to calibration uncer- 
tainties and background subtraction, was done by Hughes 
et al. (1993) for the GINGA observation of Coma. They 
found a systematic uncertainty of ±0.13 keV which, in 
fact, dominates the especially small statistical uncertainty 
(±0.09 keV) for this cluster. Arnaud et al. (1992) also found 
a typical systematic uncertainty due to background subtrac- 
tion of 0.14 keV for A2163 (T = 14.5 keV). Although this 
systematic error is expected to depend on cluster flux and 
temperature, these typical figures show that, for most of the 
clusters considered here, statistical errors still dominate the 
systematic errors. This is further confirmed if one compares 
the temperature estimates from various instruments. Gener- 
ally, good agreement between MPC, EXOSAT and GINGA 
estimates is found (David et al. 1993), as well as between 
ASCA and other high-energy instruments (Markevitch & 
Vikhlinin 1997). A direct comparison between the temper- 
ature estimate for the 7 clusters in our sample with both 
GINGA and EXOSAT measurements gives a mean ratio of 
0.97 with standard deviation of 0.06. 



3 THE L x -T RELATION AND 
IMPLICATIONS 

3.1 The Lx-T Relation 

Figure 1 shows the luminosity-temperature relation for the 
sample, with 90% confidence errors plotted for T. Statisti- 
cal errors on Lx are typically smaller than the plotted point 
size and so are not shown. The homogeneous GINGA sam- 
ple data are shown as filled circles, the remainder as open 
circles. The set of clusters define a remarkably tight rela- 
tion. The main outliers are A1060 and Triangulum, which 
are anomalously dim (or hot) in comparison to the others. 
The deviation of Triangulum is significantly reduced if the 
ROSAT luminosity is employed in place of the EXOSAT 
value. As noted above, this difference is likely due to its 
large angular extent. The peculiarity of A1060 was noted by 
Loewenstein & Mushotzky (1996) in a direct comparison to 
AWM7, a cluster of similar temperature which is also in our 
sample. Girardi et al. (1997) also pointed out the departure 
of this cluster from the general Lx — T relation. They noted 
that its temperature is high compared to its velocity disper- 
sion and suggested some anomalies in the dynamical state 
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Figure 1. The luminosity — temperature relation for the 24 clus- 
ters used in this analysis. Filled circles are the homogeneous 
GINGA subset while open circles denote clusters with ASCA 
temperatures and EXOSAT, MPC, or ASCA luminosities. Er- 
ror bars plotted on T are 90% confidence limits. The solid line 
is the best-fit relation; the dashed lines are offset by an amount 
equal to the rms deviation in log(Ljf) about the mean relation. 
The upper insert shows a histogram of the luminosity residuals 
(in dex), which makes apparent the two principal outliers, A1060 
and Triangulum. 

of the gas. However they assumed a temperature of 3.9 keV 
from David et al (1993), which is significantly higher than 
the more precise ASCA value adopted here. From the veloc- 
ity dispersion, a = 633km/s, determined by Girardi et al. 
(1997) and the best fit a — T relation established by Girardi 
et al. (1996), we actually expect a temperature of 2.8 keV, in 
reasonable agreement with the ASCA value of 3.1 ±0.2 keV. 
The deviation of A1060 is only sligtly decreased if this value 
is used instead of the ASCA value. Either both the galaxies 
and the gas are in a very particular dynamical state or the 
gas structure or content in this cluster is abnormal. 

The data are well fit by a power law \og(Lx) = (2.88 ± 
0.15)log(T/6 keV) ± (45.06 ± 0.03). The central values and 
uncertainties in the fit parameters are derived from sim- 
ple least squares fits (in log space) to a large number of 
bootstrap resamplings employing the estimated T measure- 
ment errors (quoted errors in slope and intercept are la). 
Note that since bootstrap procedure is used to estimate er- 
rors the exact fit method is of secondary importance. The 
slope of 2.88 ±0.15 is consistent with that found by Ebeling 
(1993) from RASS luminosities (2.73 ± 0.39) but is signif- 
icantly smaller than the value of 3.36 found by David et 
al. (1993) using mostly MPC and EXOSAT data. This is a 
direct consequence of the underestimate of luminosities for 
nearby clusters by these last two experiments (see above); 
the observed clusters with low luminosities being on average 
at smaller distance than more luminous clusters which can 
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be observed up to larger distances. Our slope is also con- 
sistent at better than the 2a level with the recent value of 
2.64 ± 0.16 determined by Markevitch (1998) from analysis 
of ASCA temperatures and ROSAT luminosities of nearby 
clusters. Allen & Fabian (1998) find a slope of 2.9±0.3 in the 
non cooling flow portion of their combined ASCA/ROSAT 
sample. 

The raw scatter in decimal log(Lx) about the 
luminosity-temperature relation is 0.14, and some of the 
variance arises from measurement errors in temperature. We 
estimate this contribution to be 0.06, and thus calculate the 
intrinsic scatter in log(Lx) at fixed temperature to be 0.13, 
equivalent to a fractional deviation ((SLx / Lx) 2 ) 1 ^ 2 = 0.30. 
For the homogeneous GINGA subset of 18 clusters, the frac- 
tional deviation reduces substantially, to 0.17. 

Considering A1060 as an outlier, the fact that there is 
no significant trend toward larger scatter in Lx at low tem- 
peratures implies that cooling flows do not contribute sub- 
stantially to the total luminosity of these systems. Our se- 
lection criterion of a fixed cooling flow mass deposition rate 
is roughly equivalent to an absolute limit in core luminosity 
excess attributed to the cooling flow region. This absolute 
limit represents a larger fraction of the total luminosity from 
smaller systems, leading to the concern that cooling flows of 
variable strength (but under our 100 Mq yr _1 limit) might 
induce larger variance at low temperatures if the cooling re- 
gion dominated the luminosity. This effect is not apparent 
in this data set of modest number. 



3.2 Limits on Gas Fraction Variations 

The small scatter in Figure 1 supports the notion that at 
least this sub-population of weak cooling flow clusters is 
quite regular in its structural properties. In particular, the 
data appear at odds with the idea of large gas fraction vari- 
ations within the virial regions of clusters, an idea hinted at 
by the comparison of A1060 to AWM7 by Loewenstein & 
Mushotzky (1996). We now address this issue using a spher- 
ically symmetric model as a framework for calculation. The 
purpose of the exercise is to provide a context for connect- 
ing the scatter in gas fraction at fixed temperature with the 
observed variance in X-ray luminosity. The exercise is non- 
trivial because luminosity scatter can be due to changes in 
gas density structure as well as overall gas fraction varia- 
tions. 

Consider a cluster with smooth, spherically symmetric 
distributions of gas and total mass about its center. Let 
M 0c be the total mass contained in a sphere (of radius rs c 
about the center) which encompasses a mean density 8 c p c , 
where p c = 3 H 2 /8nG is the critical density of the universe 
within which the cluster is embedded. The bolometric X-ray 
luminosity of a cluster is 



d 3 r P 2 gas {r) A(T(r)) 



(1) 



where A(T(r)) is an appropriately normalized emissivity de- 
pendent only on temperature. In practice, the integral ex- 
tends not to infinity, but to a finite radius denoting the edge 
of the relaxed, virialized portion cluster. Spherical collapse 



models and numerical simulations suggest that this edge is 
located at density contrast <5 C ~200, or radius r2oo- 

Writing the gas density in terms of the natural radial 
variable y = r/r 0c , we introduce the structure function h(y) 



Pgas(yr Sc ) = fgas S c p c h(y), 



(2) 



where the explicit use of the gas fraction within density con- 
trast 8 C 



fgas 



M ga3 (< r Sc 



(3) 



sets the normalization of the structure function h(y) via the 
condition 3 J Q dyy 2 h(y) — 1 . For an ensemble of clusters, 
we imagine a set of gas fraction values {/ 9 „} and structure 
functions {h(y)}. 

The luminosity in equation (hi) with this notation be- 



Lx = fgas Q(5 C ,T) 

where 



Q(Sc,T) 



d 3 y h 2 (y) [k{T{y))p 2 5 2 r 5a 



(4) 
(5) 



is largely determined by the structure function h(y), with 
weaker dependence on the form of the temperature distri- 
bution through T(y). 

The variance in luminosity at fixed temperature has 
contributions from gas fraction variations and structural 
variations 

Si = 43} + 4S f 6 Q + 6%, (6) 
where the S's are fractional variations at fixed temperature 

( SLx \ r _ ( Sfga 



fgas 



(7) 



The magnitude of intracluster gas variations 8j can be 
determined from the observed scatter in luminosity 8l 

S 2 f = ~ 5 2 Q ) - 5 f 8 Q (8) 

with an additional assumption about how changes in gas 
fraction and cluster structure are correlated. We assume here 
that there are no correlations between changes in gas content 
and variations in internal structure, meaning 



(8 f S Q ) = 



(9) 



with the angle brackets denoting an ensemble average. 

With this assumption, the fact that 8q is positive defi- 
nite leads to an upper limit on gas fraction variations from 
equation (0) 



(s}) 1/2 < \ {sl) 1/2 



(10) 



(11) 



and the resultant values from Figure 1 are 

/j)2o/2 < ( 0.15 full sample, 
{ f) ~ \ 0.09 GINGA subsample. 

Enlarging these limits is possible if an anti-correlation 
exists between Sf and Sq. Such an anti-correlation would 
require gas poor clusters to be more centrally concentrated 
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than average, and vice-versa. In this way, the change in lu- 
minosity due to a lower (higher) gas mass would be compen- 
sated by an increased (decreased) central density. We show 
below that the observational data do not support such an 
anti-correlation. In addition, one mechanism capable of gen- 
crating gas loss from clusters has the opposite effect. Heating 
due to winds from early-type galaxies leads to a less cen- 
trally concentrated ICM structure and slightly depressed gas 
fractions within the virial radius (Metzler & Evrard 1997). 
This implies (S/Sq) > 0, a result which serves to strengthen 
the limit imposed by equation (^). 

We conclude that variations in gas fraction within the 
virial regions (overdensities 5 C ~ 200) of weak, cooling 
flow clusters are likely to be quite small. Fractional devi- 
ations at fixed temperature 5/ may have standard deviation 
smaller than 10%, particularly if galactic wind feedback, or 
some other, non-gravitational, entropy generating mecha- 
nism, has produced correlated variations in gas fraction and 
density structure. 

The limit on gas fraction variations applies to clusters 
of similar temperature; there remains the possibility that 
the mean cluster gas fraction varies, perhaps considerably, 
over the range of temperatures probed in our sample. 



4 TEMPERATURE DEPENDENCE OF ICM 
STRUCTURE 

In this section, we examine our sample for evidence of vari- 
ation in the mean cluster ICM structure and virial gas frac- 
tion with cluster temperature. We also discuss the origin of 
the Lx~T relation slope. 

4.1 The slope of the Lx~T relation 

Suppose that cluster density profiles /ir(y) and gas fractions 
fgas(T) are, in the mean, temperature dependent. If we fur- 
ther suppose that clusters are internally isothermal, then 
equation (^) can be rewritten in a form which makes clear 
the dependence on the shape function hr(y) 

Lx{T) = f 2 gas (T) [M Sc (T)p c S c A(T)] Q(T) (12) 

with Q(T) — (3/47r) J d 3 yh^(y) is the dimensionless emis- 
sion measure. Note that Q(T) is equal to {p 2 gas ) / '{p ga s} 2 
with the angle brackets denoting the average over the clus- 
ter atmosphere. It is thus a structure factor which depends 
solely on the gas density shape and characterises the con- 
centration of the gas distribution. 

The traditional approach is to employ the following set 
of additional assumptions : (i) pure bremsstrahlung emis- 
sion (A(T)(xT 1/2 ), (ii) virial equilibrium (M Sc ocT 3/2 ), (iii) 
structurally identical clusters (Q(T) =Ci), and (iv) constant 
gas fraction (f ga s(T) = C-i). This leads to an expected scal- 
ing relation between luminosity and temperature of slope 
two 

L X {T) oc T 2 . (13) 

The fact that the slope of the observed relation is signif- 
icantly steeper than two implies that one or more of these 



assumptions does not hold in the real cluster population. 
For example, if we retain bremsstrahlung emission and the 
condition of virial equilibrium, then the observed slope of 
the Lx~T relation leads to a constraint on structure factor 
and gas fraction dependence 

f gas 2 (T)Q(T) cx r°- 88±0 ' 15 . (14) 

The assumption of pure bremsstrahlung is quite accu- 
rate above about 2 keV, but at lower temperatures, line 
emission becomes significant unless the metallicity is very 
low. Pure bremsstrahlung emission itself does not exactly 
scales as T 1 ^ 2 due to the additional Gaunt factor. However, 
this factors varies by less than 4% in the 1-10 keV temper- 
ature range. On the other hand, the line emission reaches 
about 18% of the total emission at 2 keV for a typical abun- 
dance of 1/3 the solar value. This extra line emission tends 
to flatten the slope of the Lx~T relation, as low T clusters 
are boosted in luminosity more than their high T counter- 
parts. This effect is amplified because low T clusters tend 
to have higher abundances. Given the measured abundances 
we computed A(T) for our sample (the correct value is used 
in all the following). We found that the departure from the 
T 1,/2 law is actually small: the best fit slope of the power 
law is 0.36 instead of 0.5. 

The assumption of virial equilibrium is supported by 
cosmological, gasdynamic simulations (Bryan & Norman 
1998 and references therein) . The relation holds particularly 
well for S c ~ 200 — 2500, with 200 the fiducial choice of virial 
radius and 2500 approaching the onset of the cluster core. 
We use the virial theorem as one method for estimating total 
cluster masses below. 

There remains the issue of structural regularity. Recent 
numerical investigation of dark matter halos evolved within 
hierarchical clustering cosmogonies indicates that low mass 
halos are, on average, more centrally concentrated than high 
mass halos (Navarro et al. 1997). Although the magnitude 
of the change is likely to be small within the mass range 
probed by the X-ray clusters in our sample, the qualitative 
direction is again to flatten the slope of the Lx~T relation, 
if one assumes that the gas density follows that of the dark 
matter. 

The assumption that the gas traces the dark matter is 
reasonable if shock heating through gravitational infall is 
the only entropy generating mechanism, but this assump- 
tion breaks down if extra energy sources or sinks are avail- 
able to the gas. In particular, if galactic winds play an im- 
portant role in heating the ICM, then the extra entropy 
dumped into the gas can change its density structure dra- 
matically (Metzler & Evrard 1994) . The net result of winds 
is to inflate the gas distribution, decreasing the Q factor 
and causing diminished gas densities and, correspondingly, 
lowered gas fractions, at radii interior to the virial radius. 
Because supernova-driven winds impart roughly a fixed spe- 
cific energy into their surroundings, their impact is felt more 
strongly in shallower potential wells. This steepens and low- 
ers the intercept of the Lx~T relation compared to models 
without winds. 

The magnitude of this effect will depend on details of 
the assumed wind model, but plausible models exist. The 
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Figure 2. Values of the gas outer slope parameter (3 (bottom) 
and structure factor Q (top), derived from j3— model fits to X-ray 
images plotted against temperature for 22 of the 24 members of 
our sample. Point styles are the same as the previous figures. 
Error bars available in the literature are shown. 



numerical simulations of Metzler (1994) produce a slope of 
2.96±0.05 in the bolometric Lx~T relation. Cavaliere, Menci 
& Tozzi (1997) create a semi-analytic model which predicts 
a progressive steepening of the Lx~T slope toward low tem- 
perature clusters, consistent with the data of Ponman et 
al. (1996). Finally, the constant central entropy model of 
Evrard & Henry (1991) predicts a slope of 2.75 which is 
consistent, within the errors, with our observational sample 
behaviour. 



4.2 Gas density profiles from X— ray imaging 

Azimuthally averaged profiles of cluster X-ray surface 
brightness images are generally well fit by the standard, 
/3-model form 



I(r) = Jo 1 + (r/r c ) 



21 -3/3+1/2 



(15) 



derived by assuming an isothermal gas with density distri- 
bution p gas (r) = po(l + (r/r c ) 2 )~ 3l3/2 . Fits to this model 
exist for all but two (A754 and A370) of the members of our 
sample. Table 1 lists values from the literature of the best-fit 
core radius r c and outer slope parameter /3. Errors on these 
parameters are also given in the table when available. 



Figure |2j displays the best fit slope parameters against 
temperature. As noted previously (David et al. 1990; White 
1991; Mohr & Evrard 1997), there is a clear trend toward 
decreasing (3 toward lower T. This trend is not likely to be 
an artifact of the fitting procedure, as Mohr & Evrard (1997) 
see a similar trend with f3 defined in a non-parametric, non- 
azimuthally averaged fashion. 

The heterogeneous nature and lack of error estimates for 
nearly a third of these data make it difficult to determine the 
formal statistical significance of this effect. As a reasonable 
estimate, we calculate the bi-weight modes (Beers, Flynn & 
Gebhardt 1990) for the subsamples of 8 clusters with T < 
4 keV and 14 clusters with T > 4 keV. We applied a 10% 
fractional error to each data point — slightly larger than 
the rms fractional error of 7.8% for the 15 quoted error 
values — and determined the expected subsample means 
and standard error from bootstrap analysis. The values of 
(P)io = 0.51 ±0.03 and <J3) hi = 0.66 ± 0.03 for the low- and 
high-T subsamples differ at the 3.3a level. We suspect that 
future analysis with a more regular dataset will confirm the 
statistical significance of this trend. 

The implication is that the ICM density structure of 
clusters is temperature dependent, and the sense of the ef- 
fect is such that OCT) is an increasing function of temper- 
ature (see Figure El). Models of cluster formation incorpo- 
rating galactic winds exhibit a trend of versus T similar 
to that seen in the observational data (Metzler & Evrard 
1997). The trend of Q(T) with temperature is in the right 
direction to steepen the Lx~T relation, but whether it is 
entirely responsible for this effect requires investigation of 
the trend of mean gas fraction with temperature. 

4.3 Intracluster Gas Fraction Values 

The determination of absolute values of the gas fraction re- 
quires a method to estimate the total cluster mass. We show 
here that two different methods for total mass estimation 
lead to different conclusions about the behavior of the gas 
fraction in clusters. 



4-3.1 Mass determination methods 

Estimating cluster total masses has traditionally been done 
with the "/3-model" (BM) approach (Cavaliere & Fusco- 
Femiano 1976), for which the total mass is 



Mp(r) 



1.13 x 10 15 /3 



T 



r_ (r/r e ) 2 

10 keV Mpc 1 + (r/r c )' 



M (16) 



This model assumes a non-rotating, isothermal gas distri- 
bution in hydrostatic equilibrium, with gas density profile 
of the form assumed to derive the surface brightness image, 
equation ([jl]). This leads to a binding mass which increases 
linearly with radius outside the core. Estimating the mass 
at some density threshold S c is readily done by calculating 
the mean interior density as a function of radius. Outside 
the core region, this radius scales as 



-1/2 



T 



10 keV 



1/2 



Mpc (17) 
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An alternative method is to employ the virial theorem 
(VT) at fixed density contrast, which leads to 



6- 500 



6 = 200 



M VT (rs c ) = Mio(<5 c ) 



T 



10 keV 



3/2 



(18) 



Calibration of this relation with numerical experiments by 
Evrard, Metzler & Navarro (1996, hereafter EMN) shows 
that the intercept Mio(5 c ) is remarkably insensitive to as- 
sumptions regarding the background cosmology and to the 
effects of galactic winds. The experiments show that the VT 
method has intrinsically smaller variance compared to the 
BM approach, because the imaging information introduces 
an additional source of noise. A recent intercomparison of 
12 cosmological gas dynamic techniques lends credence to 
these calibrations. Frenk et al. (1998) find agreement within 
4% for the total cluster mass and 6% for the mass weighted 
cluster temperature among codes varying substantially in 
algorithmic detail, resolving power and numerical parame- 
terization. 

At density contrast S c = 500, EMN find M w (500) = 2.2 x 
10 15 Mq, with a corresponding physical size scaling as 
r 5 oo(T) = 2.5(7710 keV) 1/2 hgj Mpc. At <5 C = 200, the mass 
intercept is Mi (200) = 2.9 x 10 15 M© and the physical 
size is r2oo(r) ^3.7/igQ 1 Mpc, roughly 50% larger than r 5 oo- 
All values are quoted for the current epoch, and both scale 
with redshift as (1 + z) _3//2 . Though typically quite small, 
this redshift scaling is taken into account in our analysis. 
We examine gas fractions at these two characteristic density 
contrasts below. For the majority of imaging data, rsoo is 
comparable to or smaller than the radius used to fit the 
gas model parameters while values at r2oo require a modest 
extrapolation outside the X-ray imaged region in nearly all 
cases. 

To determine the gas mass, we employ values of j3 and r c 
from Table 1, and determine the central gas density by nor- 
malizing to the bolometric luminosity, assuming the latter 
is due to emission within r^oo- Though not exact, the typi- 
cal error introduced by the use of V2oo as a cutoff is small, 
since for most clusters the contribution to the luminosity 
outside this radius is negligible. As a check, we normalized 
the gas mass using rsoo — an extreme choice since many 
cluster images extend beyond this radius and since rsoo is 
nearly a factor 2 smaller than r2oo — and the fractional in- 
crease in gas density was typically ~5%, with a maximum 
fractional increase of 24% (absolute increase ~ 3%) for low j3 
clusters. As detailed below, our procedure produces a pop- 
ulation mean gas fraction at rsoo for hot clusters consistent 
with the determination of Evrard (1997) from the observa- 
tional samples of White & Fabian (1995) and David et al. 
(1995). 



4-3.2 Gas fraction estimates 

Gas fractions and total masses derived from the VT and BM 
methods at 5 C = 500 and 200 are shown in Figure ^[ There 
are trends apparent which reflect the different behavior of 
the two mass estimate methods to the input data, partic- 
ularly the sensitivity to the slope parameter (3. Comparing 
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Figure 3. Gas fraction versus total mass at inferred density 
contrasts 5 C = 500 and 200 derived using the VT (upper) and 
BM (lower) approaches discussed in the text. Symbol types are 
identical to those in Figure 1. 



equations (^6|) and ( |l8| ) at rsoo, and ignoring (r c /rsoo) 2 , one 
finds 



Mv T (r 5 oo) V0.79 



3/2 



(19) 



From Figure |2j, all but one of the values of f3 in our sample 
lie below 0.79. The /3-model estimates are therefore consis- 
tently lower than those of the virial theorem, and this is 
clearly reflected in the overall horizontal shift of the data 
points between the upper and lower rows. Furthermore, the 
systematic trend of f3 with temperature means the discrep- 
ancy is larger at smaller temperatures/masses. This effect 
also causes the gas fractions at constant 5 C in the BM to 
be consistently higher than those of the VT, although the 
difference is small at S — 200. 

In some of the panels, a trend of mean gas fraction with 
cluster mass or temperature is evident to the eye. The lack 
of rigorous, homogeneous, statistical error estimates for the 
r c and f3 values in our sample makes it difficult to make pre- 
cise statements regarding the significance of these trends. We 
proceed with an approximate treatment in which we assign 
a la fractional uncertainty of 30% to each of the gas frac- 
tion estimates. This approach allows us to explore trends in 
the data and provides an estimate of their statistical signif- 
icance under an assumed error budget. The 30% error may 
be generous given that T, r c and /3 are typically determined 
to ~ 5% accuracy, and the fact that gas mass errors are 
smaller than errors in /3 and r c alone, because these param- 
eters are tightly correlated. However, the "cosmic variance" 
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Figure 4. Bi-wcight mean gas fractions at density contrasts of 
500 and 200 for the VT and BM methods. Data for the full sample 
of 22 clusters are shown as filled circles, the low T subsample of 8 
clusters as triangles, and the high T subsample of 14 clusters as 
filled squares. Error bars are 90% confidence limits (1.65cr) deter- 
mined from bootstrap resampling and assuming a fixed fractional 
error of 30% on individual cluster gas fraction estimates. 

Table 2. Bi-weight Mean Gas Fractions (%, h~ 3/2 ) 

5 C = 500 5 C = 200 

Sample VT BM VT BM 

All 14.0 ±1.3 18.3 ±1.5 18.1 ±1.4 21.0 ±1.8 

T < 4 keV 10.2 ±1.5 15.3 ± 2.4 14.5 ±2.1 19.9 ± 3.2 
T>4keV 16.2 ± 1.5 20.0 ± 1.9 20.1 ±1.6 21.6 ± 2.2 



discrepancy between A1060 and the mean of the cool sub- 
sample is smaller under the VT approach. 

The value of 16.2 ± 1.5% for the hot subsample in 
the VT case is consistent with the 95% confidence region 
fgasirsoo) = 17.0 ± 1.7% ( /150 = 1) determined by Evrard 
(1997) from the combined samples of White & Fabian (1995) 
and David, Jones & Forman (1995). Only 3 of the 26 clus- 
ters in the combined study have T < 4 keV (and these have 
fairly uncertain gas fractions which make negligible contri- 
bution to the sample mean), so comparison with our high-T 
subsample is reasonable. The agreement between these two 
studies is not entirely trivial. Although there is nearly 50% 
overlap with our subsample — 6 of our 14 are in common 
with White & Fabian — we employ here different temper- 
ature estimates in many cases (though all are consistent 
within the quoted errors) and we normalize the gas mass 
differently. 

A less obvious trend in the data is the larger variance 
in gas fractions in the BM case compared to the VT. This 
is more apparent in the T-dependent subsamples, since the 
trend of mean gas fraction with temperature mimics vari- 
ance in the full VT sample. The lower scatter in the VT 
method confirms theoretical predictions that the j3 param- 
eter acts as an extra source of noise in the BM mass esti- 
mates (EMN). An increase in scatter is also seen at larger 
radii/lower densities, consistent with the numerical model 
predictions and interpreted simply as reflecting the longer 
relaxation timescales in the outer regions of clusters. 



in the binding mass estimates is believed to be ~ 15 — 30% 
at these density contrasts (EMN), and this represents the 
limiting accuracy of any dataset. 

Bi-weight mean gas fractions determined from a boot- 
strap analysis assuming fixed 30% fractional errors in each 
measurement are listed in Table 2 and summarized in Fig- 
ure ^. Data for the whole sample as well as cool and hot 
subsamples are provided. Uncertainties quoted in the Table 
are standard errors of the mean, while 90% confidence errors 
(assuming Gaussian statistics) are plotted in the Figure. 

The trend of rising gas fraction with increasing radius 
(lower 8 C ) is apparent with both mass estimation meth- 
ods, and reflects the more extended nature of the gas den- 
sity compared to the assumed or inferred dark matter pro- 
files. Gas fraction values are consistently higher in the BM 
approach, due to the smaller, inferred total masses noted 
above. As a result, significantly different sample means can 
be quoted depending on the mass method and extent of ra- 
dial coverage of the dataset. Compare, for example, f gas = 
0.140T0.013 for the VT at <5 C = 500 with f gaa = 0.210±0.018 
for the BM at <5 C = 200. 

The range becomes wider when one considers subsets 
of the data defined by temperature. All the panels in Fig- 
ure (Q) display a trend of decreasing gas fraction at lower 
T, but only the VT method at 5 C — 500 is deemed statis- 
tically significantly (2. 8a) under our 30% error hypothesis. 
Due to the robustness of the bi-weight mean estimator, this 
conclusion is valid with or without the inclusion of the prin- 
cipal outlier A1060. ft is perhaps also worth noting that the 



4.4 Implications 

There is strong evidence from the slope of the Lx~T re- 
lation and from the behavior of the image fall-off parame- 
ter j3 that cluster structure varies systematically with tem- 
perature. Whether the virial gas fraction varies with clus- 
ter temperature is a question whose answer is both model 
and scale dependent. The model dependency we highlight 
is the choice of total mass estimation method, while the 
chosen density contrast sets the scale. For example, the 
/3-model applied at S c = 200 produces similar mean values in 
high and low temperature subsamples — 0.216 ± 0.022 and 
0.199 ± 0.032, respectively — whereas the virial theorem at 
S c = 500 yields significantly different values in these subsam- 
ples — 0.162 ±0.015 against 0.102 ±0.015, respectively (see 
Table 2 and Figure ^, all values scale as hso~ 3 ^ 2 )- 

That the implications of the observed luminosity- 
temperature relation are model dependent can be seen by 
considering the behavior of terms with temperature depen- 
dence contained in the denning relation, equation ( [l2] ) . The 
structure factor Q(T) appears to increase with temperature, 
but this effect is not solely responsible for the steepening of 
the Lx~T relation. The dependence with T of the other 
factors depends on the total mass estimation method. In 
the VT approach, while Ms (T) scales classically as J l3 / 2 ] 
fgas{T) increases with T. In the BM case, while f gas is found 
to be nearly constant, the slope of the M$ c (T) — kT rela- 
tion is steeper than 3/2, due to the extra /3 factor and its 
increase with T (equation (^)). Both cases result in a steep- 
ening of the fg a3 (T) Ms c (T) term to a slope greater than 
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the canonical 3/2. The temperature dependence of this term 
contributes about equally with Q(T) to the steepening of the 
Lx-T relation. 

Allan & Fabian (1998) take a different approach in anal- 
ysis of a cluster dataset combining ASCA spectroscopy with 
ROSAT imaging. They fit ASCA data to a multi-component 
spectral model which includes emission from a central cool- 
ing flow (Johnstone et al. 1992), absorption from galactic hy- 
drogen, intrinsic hydrogen absorption, and a superimposed 
background plasma at temperature T. Compared to simpler 
fits without cooling flow and intrinsic absorption, the best fit 
temperatures in their complex fits are higher, and the best 
fit slope of the Lx~T relation is shallower, 2.33 ±0.43. Allen 
& Fabian emphasize that this slope is consistent with the 
straightforward theoretical expectation of 2; however, note 
that the range of slopes allowed in our fit does overlap the 
error range quoted above (derived from a change 8\ 2 = 2.71 
in their \ 2 analysis). 

The approximate nature of our error analysis limits the 
conclusions which can be drawn regarding gas fraction vari- 
ations. A definitive analysis awaits uniform treatment of a 
large, homogeneous dataset. There is some supporting ev- 
idence for deficient gas fractions in poor groups (Hickson, 
1997 and references therein), but there also appears to be 
larger scatter at low temperatures (e.g., Ponman et al. 1996) 
and observational uncertainties on the ICM mass remain 
large in these systems. In particular the separation of the 
ICM emission from the X-ray emission linked to individual 
galaxies is ambiguous and difficult with present instruments 
(Mahdavi et al. 1997, Mulchey and Zabludoff 1998). 



5 CONCLUDING DISCUSSION 

The narrow scatter in the Lx~T relation of weak cooling 
flow clusters (M < 100 Mq yr _1 ) indicates that the struc- 
ture of the intracluster medium in this majority sub-class 
is very regular. Under a conservative assumption of no cor- 
relation between virial gas fraction and internal, structural 
variations {(SgSf) = 0), we place upper limits on the rms 
percentage variation in virial gas fraction of 15% for the 
entire sample of 24 clusters, and 9% for the homogeneous 
GINGA subset of 18 clusters. 

These limits are conservative because the X-ray imag- 
ing data supports the point of view that gas in low T clusters 
is more extended than that of high T clusters (Figure |^). 
Mohr & Evrard (1997) find the same result with a non- 
parametric analysis of X-ray images. They also derive limits 
on gas fraction variations similar to those found here based 
on the tightness of the X-ray size-temperature (ST) rela- 
tion. The agreement is non-trivial because the ST analysis 
is essentially local, being based on an differential isopho- 
tal area, while the Lx~T relation analysis is global, based 
an the integrated emission of the entire cluster atmosphere. 
Both analyses are consistent with idea that energy input 
from galactic winds has led to extended gas atmospheres 
in low T clusters (David, Forman & Jones 1990,1991; White 
1991; Metzler & Evrard 1997; Cavaliere et al. 1997; Ponman 
et al. 1996). 



In wind models, gas loss is correlated with structural 
extension of the gas, leading to (<5q<5/) > 0. The inferred 
variance in gas fraction derived from equation (^) could be 
substantially smaller than the limits phrased above if winds 
have played an important evolutionary role. Numerical sim- 
ulations of clusters used by Mohr & Evrard (1997) exhibit 
less than 5% fractional variation in their virial gas fractions, 
and a similar value is seen in models with galactic winds 
(Metzler & Evrard 1997). The expected "cosmic variance" 
in virial gas fractions is thus quite small. 

Comparison of a number of different simulation tech- 
niques paints a similar picture. Frenk et al. (1998) find agree- 
ment at the 5% level (la deviation) for the virial gas fraction 
of a single cluster determined by twelve independent gas dy- 
namic methods. This indicates that systematic uncertainties 
in the tested numerical simulation techniques are small for 
the simplest case of a single phase ICM. However, it is im- 
portant to note that the quoted simulations ignore galaxy 
formation and its attendant multi-phase, gas structure. In- 
troducing galaxy formation via cooling and star formation 
is likely to increase the scatter in ICM gas fraction values, 
and it will be interesting to see whether viable models can 
do so without violating the observational limits presented 
here. 

Evidence for trends in f gas with cluster temperature is 
ambiguous. A key uncertainty lies in estimation of the to- 
tal, binding mass. The /3-model approach produces masses 
which are characteristically smaller than the virial theorem 
method as currently calibrated by numerical simulations. 
For a cluster with f3 = 0.6, the masses inferred at S c = 500 
differ substantially, Mp = 0.QQMvT- This different behavior, 
coupled with the temperature dependence of (3, leads to dif- 
fering conclusions on whether the mean gas fraction varies 
with cluster temperature. 

For the same reason there is some ambiguity in the ori- 
gin of the steepening of the Lx~T relation. The fact that 
the gas is more concentrated in high-kT cluster explains at 
least half of the effect, directly through the structure factor 
in the X-ray luminosity. In the /3-model approach, the tem- 
perature dependence of the ICM structure can fully account 
for the Lx~T relation slope, via the additional dependence 
of the total mass with (3. On the other hand, in the VT 
approach, the gas mass fraction increases with T and this 
effect contributes about equally to the Lx~T relation slope 
increase. Both effects, lower gas fraction and inflated gas 
distribution in low T systems, are expected in models incor- 
porating the effects of galactic winds. 

Independent methods for estimating virial masses are 
crucial, and weak gravitational lensing provides an impor- 
tant approach. Allen (1997) provides evidence that weak 
lensing and /3-model mass estimates are consistent, but the 
present errors at radii between r$oo or r2oo are too large to 
discriminate between the BM and VT approaches. We sug- 
gest that weak lensing analysis of a small sample of mod- 
erate redshift clusters selected to have /3 < 0.6 from X-ray 
imaging could settle the score, if individual measurement 
uncertainties could be kept below a few tens of percent at 
<5 C = 500. 

Future X-ray missions which provide simultaneous 
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imaging and spectroscopic capability (AXAF, XMM, Astro- 
E) will provide improvements over the current dataset. The 
ideal analysis would be one of a flux limited sample with 
central, cooling flow regions excised (Markevitch 1998), us- 
ing a nearly homogeneous dataset of temperatures, lumi- 
nosities and surface brightness maps derived from the same 
telescope. Such data would minimize variance caused by in- 
strumental uncertainties and open the door to very high 
precision measurements of intracluster gas fractions. 

Finally, small scatter in virial gas fractions lends sup- 
port to arguments against an Einstein-deSitter universe 
based on cluster baryon fractions and primordial nucleosyn- 
thesis (White et al. 1993; Evrard 1997). The standard coun- 
terargument to this idea has been that the physics of the 
ICM is more complicated than is assumed in most analytic 
or numerical models. However, if additional physics beyond 
gravity, shock heating and galactic wind input is at work in 
cluster ICM atmospheres, then the extra processes involved 
must conspire to keep the scatter in the observed Lx~T and 
ST relations small. Since Lx and Ri (the "size" in the ST 
relation) are manifestly different measures (global vs. local) , 
it is not obvious that interesting physical mechanisms, such 
as a strongly multi-phase ICM (Gunn & Thomas 1996), can 
be added in such a way as to preserve the current agreement 
between observations and current models incorporating the 
effects of galactic winds on the ICM. 
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